Temps = [0.35,0.35,0.45,0.6,0.8,1.3,4,8.5];
TempsErr = [0.1,0.1,0.1,0.15,0.1,0.3,0.5,0.5];
%TODO: get reasonable temp error bars

Ds = [37.9,37,26.8,26.7,21.5,18.1,15.2,13.3];
DsErr = [2.4,4,1.5,1.6,1.8,1.2,0.9,0.8];
Gammas = [4.3,4,6,6.1,6.5,6.2,7.2,11.8];
GammasErr = [0.4,0.4,0.9,0.5,0.8,0.5,0.9,2.3];
Visc = [6.4,8.2,3.6,2.5,2.9,2.7,2.4,8.5];
ViscErr = [0.9,1.4,0.6,0.9,0.9,0.4,1,2.2];

figure;
subplot(2,2,1)
errorbar(Temps,Ds,DsErr,DsErr,TempsErr,TempsErr);
grid on;
xlabel('T (t)');
ylabel('D (latt^2/ms)');
ax = gca;
ax.YLim(1) = 0;

subplot(2,2,2)
DInvErr = DsErr./Ds.^2;
errorbar(Temps,1./Ds,DInvErr,DInvErr,TempsErr,TempsErr);
xlabel('T (t)')
ylabel('1/D ms/latt^2');
grid on;
ax = gca;
ax.YLim(1) = 0;

subplot(2,2,3)
errorbar(Temps,Gammas,GammasErr,GammasErr,TempsErr,TempsErr);
grid on;
xlabel('T (t)');
ylabel('\Gamma (KHz)');
ax = gca;
ax.YLim(1) = 0;

subplot(2,2,4)
errorbar(Temps,Visc,ViscErr,ViscErr,TempsErr,TempsErr);
grid on;
xlabel('T (t)');
ylabel('\nu (latt^2/ms)');
ax = gca;
ax.YLim(1) = 0;

suptitle('Model Params Vs. T')


